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Abstract 

The possibility of separating the 7r'^,r/, and meson background from the signal photons 
produced directly in pp collisions is analyzed. The rejection factors for two pseudorapidity regions 
0.4 < r/ < 1.0 and 1.6 < < 2.0 and six Et values, 20, 40, 60, 80, 100 and 200 GeV, are calculated 
for a case when only the calorimeter information is used. The cases of r/, a; and K'^ mesons 
decaying through the neutral and charged channels are considered separately. 



1. Introduction. 



This work is a continuation of our previous publications on the study of "7 + 
jet" events at LHC energies [|l]] - [||]. In those papers it was shown that a process of 
direct photon production (with Pp > 40 GeV/c) in association with one jet, caused 
mainly by the Compton-like qg ^ ^ + q and annihilation qq ^ j + g subprocesses, 
has a considerable background due to other QCD processes that contain high Pt final 
state photons originating from decays of vr", t], Kg and oj mesons (see The former 
subprocess, as was pointed out in [Q], can be used at the LHC for studying the gluon 
density in a proton in the reaction of inclusive direct photon production and in the 
"7 + jet" events ([||], [p^], [11]). As for the decay channels of the abovementioned 
mesons (see Table 1 with the PDG data [13] of branching ratios), we shall consider 



first the neutral decay channels (vr*' — > 27; —>■ 27:^; rj — > Stt^ and 27) as the next 
most important background source (especially with increasing signal photon energy) 
after the processes with a hard radiation of photons from quarks, i.e. bremsstrahlung 
photons (see [^). 

Table 1: Decay modes of vr'^, rj, and uj mesons. 
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We shall study a question of what level of accuracy can be achieved in suppression 
of the contribution from these neutral decay channels if we confine ourselves only to 
the information from the CMS electromagnetic calorimeter (ECAL). Possible use of 
preshower detector information is not considered here as it was a subject of another 



publication [|12|]. 

Another group of decay channels that contain charged pions in the final state 
is less difficult to be suppressed. As will be shown below, their contribution to the 
background can be discriminated with a good efficiency on the basis of the hadronic 
calorimeter (HCAL) data. 

The results presented here are obtained from the simulation with the GEANT 
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based package CMSIM f| We carried out a few simulation runs including 

(1) four particle types: single photons 7 and vr*^, r/, mesons; 

(2) six Et values for each of these particles: 20, 40, 60, 80, 100 and 200 GeV; 

(3) two pseudorapidity regions: 0.4 < 77 < 1.0 (Barrel) and l.6<r]< 2.0 (Endcap). 
About 4000 - 5000 single-particle events were generated for each Et, rj interval and 
for each type of particle. 



2. Neutral decay channels. 

To separate the single photon events in the ECAL from the events with photons pro- 
duced in the neutral decay channels of vr'', r/, and lo mesons, two variables which 
characterize the spatial distribution of Et deposition in both cases were used: 

1. Et deposited inside the most energetic crystal cell in the ECAL Et'j^f^^', 

2. The quantity D^j that characterizes spatial Et distribution inside the 5x5 crystal 
window around the most energetic crystal cell (see below). 

The showers produced in ECAL crystal cells by the single photon are concen- 
trated in a 5x5 array of crystals centered on the crystal with the maximum signal 
(see [16]). Another picture is expected for the multiphoton final state arising from the 
meson decay through the intermediate tt^ states that would cause different spatially 
separated centers of final state photon production. In this case the spatial distribu- 
tion of Et deposited in such a shower is also more likely to be different from the one 
produced by a single photon. 

Let us define the coordinates of the ECAL shower center of gravity (gc) in the 
T] — (f) space according to the formula: 





V9c=\^ mEni \l^Et^\; </),c = 1^ <I^.Et^\/\ 1^ Et^ (1) 



The sum in formula (|I|) runs over ECAL crystal cells forming the N x N window that 
contain a shower with the most energetic cell in the center. Here i?t* is Et deposited 
inside the ith crystal cell belonging to the shower. 

The distributions of Et deposited in ECAL by the single photon, 77 and vr'^ 
mesons with initial Et = 20, 40, 60 and 100 GeV are shown in Fig. [I] as a func- 
tion of the distance counted from the center of gravity Rgc- One can see from Fig. |l| 
that about 96% of the total photon Et deposited in ECAL, i.e. Et-tot, fits into the 
radius Rgc = 0.02 and practically all its energy is contained inside Rgc = 0.045. 
From the same pictures we can conclude that for all mesons the total Et of showers 

' The results in Sections 2 and 3 of this paper are obtained with CMSIM versions 111 and 116 respec- 
tively. 

^(with acell size of 0.0175x0.0175) 
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is mainly contained inside the 5x5 crystal cell window having the size of one CMS 
HCAL tower. So, in formula ([I]) we can put with a good accuracy N = 5. We also 
observe from this figure that the difference in spectra for single photons and pions 
decreases with growing Et and practically disappears at Et = 100 GeV, while the 
difference of the distributions caused by rj mesons from those for single photons can 
still be seen at Et = 60 and 100 GeV. 
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Figure 1: The normalized dEt/dRgc/ Et^tot distribution over Rgc of the showers 
produced in ECAL by vr*^, 77 mesons and single photons 7 having Et = 20, 40, 60 and 
100 GeV. 

For our further needs we introduce another useful variable. Let us consider the 
distance of the ith cell from the center of gravity of the electromagnetic shower in 
ECAL, i.e. 

m-<t^9c? + {m-Vac?Y'\ (2) 



where (r^j, 0j) are coordinates of the center of the ECAL crystal cell. The distances 
Tj of each crystal cell can be used to introduce a new useful quantity (that effectively 
takes into account the contribution of energetic cells far from the center of gravity): 

25 



A, 




(3) 
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For illustration, we present the results of CMSIM simulation with E{^^^^^ = 
40 GeV (mes = 71*^,7/) in Figs. ^ and|3[ These figures contain the normalized distri- 
butions of -Etmax ^rid Djj^ that characterize the showers produced by the products of ij 
meson (Fig. ^ and vr*^ meson (Fig. ^ decays. The analogous Et^^^^ and distribu- 
tions in the showers produced by a single photon are also shown there for comparison. 
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Figure 2: The normalized distributions of the number of events over and Et^^j^^ for 
the single photons (7) and for ECAL showers originating from r/ mesons (77) having 
Et = 40 GeV. 
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3: The normalized distributions of the number of events over 
single photons (7) and for vr*^ mesons (vr*^) having Et = 40 GeV . 



and Ei'^^}^ 



One can easily see that the 77 meson shower spectrum over is strongly dif- 
ferent from the one of the photon. The range of Et^^^^ values where the 77 meson 
background makes a small contribution {Et'^l^^ > 28 GeV) is also clearly seen in 
Fig. ^ From the result of the simulation we conclude that 90% of the signal pho- 
ton events at SjT/"^^^ = 40 GeV are concentrated in the regions < 0.0088 and 
> 23 GeV. 

So, the Et^l^ and Du, spectra of the i] meson can be effectively used for sepa- 
ration of the background and the signal events because the region of their overlapping 
becomes much smaller as compared with Fig. 1 . 
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We refer to the way of separation of signal and background events based on 
the differences between distributions over the 1st variable -Etmix "Et-max" 
criterion and over the 2nd variable as the ''D^w" criterion. 

The situation is much worse in the case of vr'' showers as is seen in Fig. |3| where 
the spectra of practically overlap with the spectra of the photon. 

The neutral pion rejection efficiencies R^jf (relative to direct photons) ^ ob- 
tained on the basis of the D^w criterion for different Ep/'"° are given in Table ^ as a 
function of the chosen single photon selection efficiencies (in %) eejj= 80, 85, 88 and 

90 %. The second line in this table ''E{^l'"^ = 40 GeV" corresponds to the plots in 
Fig. |. We see from this table that -Re// grows by almost two-fold as ejjj decreases 
from 90% to 80%. 

Table 2: Neutral pion rejection efficiencies -Rgjj (%) obtained from CMSIM simu- 
lation by application of D^w criterion for five Et values of single 7 and vr*^ and for 
different values of single 7 selection efficiencies r r = 80—90%. 0.4 < rp^'^ < 1.0. 
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The rejection percentages of events with vr for each of the criteria Et^max 
and D-w are given in Fig. 4 for different meson Et Q 

Analogous rejection curves obtained for the rj and Kg meson neutral channels 
(for the fixed value of single photon selection efficiency ejjj- = 90%) are presented in 
Fig. 5 (Et^max criterion) and in Fig. 6 (D^w criterion) for the interval 20 < Et < 
200 GeV. 

By comparing Fig. 5 for the rj meson and Fig. 4 we see that for the same 
^2ffi— 90%) and the same rejection criterion Et-max the vr'^ rejection efficiency 

Reff becomes less than 30% at Et = 30 GeV while the 77 meson rejection efficiency 
drops to the same level of 30% only for Et > 90 GeV. It is also seen that 



''The rejection efficiency is defined as a ratio of the number of background events discarded by the 
cut, taken for a given value of signal selection efficiency e^//' '° number of background events. 

''The rejection factors for the Endcap are found to be in agreement with those of |jl^ (without 
preshower) if one takes into account the fact that in Fig. 4 the rejection powers are averaged over the 
entire 1.6 < 77 < 2.0 range. 
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D^w criterion works a bit better than Et_max one for the Barrel region (0.4 < r] < 
1.0). 

In spite of a difficuhy of neutral pion background separation from single photon 
signal it is worth reminding that, as it was already mentioned in the vr*^ background 
events contribution (as well as the contribution from r/, lo and events), left after 
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Fig. 4: TT rejection efficiencies for Et-max and D-w criteria (e^ 
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the cuts chosen in and [0] decreases with growing faster than the "7 — brem" 
contribution. For this reason the bremsstrahlung photon background is more dangerous 
than the one from vr^ events. It was also shown in that the contribution to the 
total background from the events containing the uj meson as a candidate for the direct 
photon is less than 1 — 2%. Besides, the part of the neutral decay channels is only 8.5% 
(Table 1). That is why we have not considered in this section the rejection possibility 
of the Lu mesons decaying via neutral decay channels. 

3. Charged decay channels. 

Now let us return to second group of the mesons decay channels that have charged 
daughter particles in the final state (28% of t], 68.6% of and 88.8% of lo meson 
decays; see Table 1). 

We have studied the angle separation between the charged pions originating 
from LO and t] mesons and the neutral pion (as in a case ri,uj —>■ vr+Tr^vr") or photon 
(rj — > 7r+7r~7) produced in the same meson decay. For this aim the simulation of 1 
million pp events at ^/s = 14 TeV was carried out using PYTHIA 5.7 with the set of 
all QCD and SM subprocesses having big cross sections. The minimal Pt of the hard 
2^2 subprocess, i.e. p ™™ = CKIN{3) parameter in PYTHIA, was taken to be 
prnm _ Q^y^ ^ j^e Corresponding spectra normalized to unity are presented in 
Figs. |7| and I for ^i'^''' > 30 GeV separately for the charged pions with the maximal 
and minimal Et. One can see (Fig. ^ that the charged pion deflects from the vr*^ 
direction in the hadronic decay of to by the angle A9 = 0.4 — 0.5° and by the angle 
A(j) = 1.0 — 1.3°, on the average (the (p size of one crystal cell is about Acj) = 1°). The 
corresponding averaged deflections in the abovementioned hadronic decays of the rj 
meson are = 0.2° and Acj) = 0.5° (see Fig. ^). Thus, from the distributions shown 
m Figs. and I we may conclude that practically in all events the charged pions enter 
the same 5x5 EC AL crystal window as vr*' (7) does and they may partially deposit 
their energy in the HCAL towers behind the ECAL cells that register the 7r°(7) signal 
0. The appearance of the corresponding signal in HCAL (see below for details of the 
CMSIM simulation) may allow, in principle, rejection of the events with the hadronic 
background. 

The Et spectra of the charged pions originally produced in these channels were 
also studied using PYTHIA. They are presented in Fig. ^ for parent r], and lo 
mesons having Et = 40 GeV separately for charged pions with maximal Et'^^^ and 
minimal Et^i^ decay event. We see that the Et^^x distribution spectra for 



^Strictly speaking, these results are valid only on the PYTHIA level of simulation and may be slightly 
modified when the magnetic field effect is taken into account in full GEANT simulation with CMSIM. 
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mesons having Et = 40 GeV start at about 17 — 18 GeV for vr^ from decays and 
at about 5 — 6 GeV for tt^ produced in decays of rj and uj mesons. 

So far we were discussing the spectra obtained at the level of PYTHIA simula- 



tion without the account of detector effects. In Fig. 10 we present two spectra obtained 
after CMSIM simulation of the calorimeter response to the propagation of charged pi- 
ons in the CMS detector (the complete CMS setup was used). They correspond to the 
case of charged pions having Et = 5 GeV and pseudorapidity r] equal to 0.4 and 1.7. 
The spectra are normalized to the total number of events and describe the distribution 
of transverse energy deposited in the HCAL (''Et^dep"). They are built using 2000 
CMSIM simulated events with a single charged pion for each rj direction. 
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Figure 7: Absolute values of the difference in the 9 and cf) angles between vr^ and 
photons or vr'^ originated from the uj vr+vr^vr'^ decay. Plots la (lb) and 2a (2b) 
correspond to the spectra over the angles between and the charged pion with the 
maximal(minimal) Et'^ in the decay channel. 
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One can see from Fig. [T^ that even a single charged pion with Et = 5 GeV has 
enough deposited energy to produce a noticeable signal in the HCAL in more than 90% 
of events. In decays mentioned in Table 1 the charged mesons are produced in pairs. 
The Et spectra of the meson having the largest Et in this pair start at Et^^^ ^ 5 GeV 
as may be seen from Fig. ^ So, one can expect that at least 90% of this background 
may by rejected by measuring the HCAL signal. 

As we noted above, in reality we shall have a combined contribution of two 
charged pions to the 5x5 ECAL crystal cell window in one decay (see Figs. and ||) 
and, thus, to HCAL towers behind it. The following CMSIM simulations were carried 
out to investigate this contribution. First we considered the rj meson decay 
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channel. 
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Figure 9: Et spectra of charged pions from K^, r] and lo meson decays. Dashed and 
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Figure 10: isf deposited in the HCAL by charged pions with initial Et = 5 GeV. The 
left-hand plot corresponds to the Barrel region, the right-hand one corresponds to the 
Endcap region. 
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as a typical example. To examine in what way the charged mode of t] meson decay 
can fake the signal photon we forced rj meson to decay only through the charged 
channels (rj tt^tt^tt^, tt^tt^^) with the pseudorapidities t] = 0.4 (Barrel region) 
and 7] = 1.7 (Endcap region). Three initial t] meson Et ranges were chosen (for both 
pseudorapidity values): -Et^„.^= 40-^60, 60^80 and 80^100 GeV. About 5000-7000 
single r] meson events for each r] direction and Et interval were generated for this aim 
(again with the complete setup of the CMS detector). 

Since we want here to separate the rj meson background one the basis of the 
presence of its hadronic decay products around a photon or a neutral pion (see Table 
1), we apply here the isolation criteria formulated in [^. It turned out that in the 
abovementioned generated samples of rj meson events the number of events with the 
transverse energy inside the 5x5 ECAL crystal cell window -EjIj^^^ > 40 GeV is 
not quite sufficient for our analysis. As a good approximation to the value declared in 
[|1 (Pp > 40 GeV/c) we choose here Et]^^^^ > 35 GeV as a lower cut. 

The spectra of the hadronic transverse energy deposited in HCAL towers behind 
the 5 X 5 ECAL crystal cell window Etjj"^ are given in plots (al), (bl) and (cl) of 
Figs. 11 and for the Barrel and the Endcap, respectively, for three abovementioned 
Et ranges. In these plots Et^ff"^ is defined as a sum of the transverse energies Et 

EtW"' = J2Et' (4) 



deposited in each ith HCAL tower fired by vr showers. No events were found with Et 
deposited in the HCAL with Et'j^"" < 0.5 GeV. In the intervals Etl^^ = 40 ^ 60 and 
60 -^ 80 GeV only about 0.1 - 0.2% of the events passed the cut £^t|^^^^>35 GeV 
have Etf"" < 1 GeV. No events with Et'^'^ < 1 GeV were found for Et^^^ = 
80 100 GeV for both values t] = 0.4 and 1.7. 

Let us remind that in we have chosen the value of the isolation cone ra- 
dius around a 7'^*'^-candidate to be i?/^^; = 0.7. It is useful to find out what radius 
of the hadronic energy shower from the rj meson decay may be in reality. For this 
aim we calculate Fh{R) = Et'j^i'^iR) / Et^E^AL- The value of Etf'^iR) differs from 
^^sum (jef^j^ejj 5y including to the sum only those towers that fit into the cir- 

cle of some radius R{r], cj)) counted from the most energetic ECAL crystal cell of the 
abovedescribed 5x5 ECAL crystal cell window, i.e. 

EtW"'{R) = J2Et'. (5) 

The dependence of Fh on this radius R describes the Et saturation of the space around 
the most energetic ECAL cell in the 5x5 ECAL crystal window. This dependence 

^containing Et^^^\^ > 35 GeV 



11 



is shown in Fig. |13| by the dashed line for three ranges of the initial transverse ener- 
gies and two values of t]. The quantity reaches the saturation values in the range of 
R{r], (p) 0.3 — 0.4 that has a meaning of a real radius of the hadronic energy de- 
position in the HCAL for the single r] meson decay. This value agrees with the value 
<o« = 0-7 chosen in [§. 

Now let us see what size of an additional Et would be added to the isolation 
cone by the energy of charged pions deposited in ECAL cells surrounding the 5x5 
ECAL crystal window (containing Et'^^j^j^ > 35 GeV). Let us define the value of 
^t^E+H ^ scalar sum of Et deposited inside the calorimeter ECAL-i-HCAL cells 
around the ECAL 5x5 crystal cell window (i.e. with subtraction of Et deposited 
inside the ECAL 5x5 crystal cell window itself): 

T-i sum— _ Sr^ rp i p 5x5 /^n 
EtE+H - 2^ - EtECAL- <.0) 

i€ECAL+HCAL 



Plots (a2), (b2) and (c2) in Figs. [11| and |12| include the spectra of Et deposited 



in the ECAL-i-HCAL calorimeter cells that are beyond the 5x5 ECAL crystal cell 
window containing the most energetic crystal cell at its center, but within the radius 
R{r], (p) = 0.7 counted from the center of this cell. They are given again for three 
different ranges of Et^^^^t ^^'^ ^^'^ values of pseudorapity t] = 0.4 and 1.7. We see 



that all spectra in Figs. [11| and [12| steeply go to zero in the region of small Et values. 
It allows the background from the charged t] meson decays to be reduced by limiting 
Et deposited in the ECAL-i-HCAL cells within the radius R{'q,4)) = 0.7. One can 



see from the right-hand columns of Figs. 1 1 and M3 that there are no events with Et 



less than 2 GeVfov the entire range of Etl^^t ^0 GeV and both pseudorapidity 

values. This fact partially explains our choice of the direct photon isolation cut in [^] 
with Et'""^ < 2 GeV/c. 

Let us introduce another variable Fe+h{R) which has a meaning of the ratio of 
EtE+H {R), i-e- the value of EtE+u taken in analogy with (|) for ie R,to Et^ECAL- 

Fe+h{R) = Et%THiR)/Et%'AL- (7) 

Its dependence on the distance from the center of gravity R{r], (p) (counted from the 
most energetic ECAL crystal cell) is shown in Fig. [l^ by solid lines. This variable 
also reaches its saturation in the same range of R{r], 0) 0.3 — 0.4. The difference 
between solid and dashed curves defines the percentage of transverse energy deposited 
by electromagnetic showers produced by charged pions from the r] meson decay in 
ECAL cells around the 5x5 ECAL crystal window containing the 7-candidate with 
EtECAL ^ 35 GeV. This difference does not exceed 14 - 20% in the Barrel and 
7 — 9% in the Endcap. Its dependence on the distance in the region of small R{r],(j)) 
is evident from these pictures. 
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Multiplying the ratios Fh{R) and Fe+h{R) in the region where they reach 
their saturation by Et^j^^j^^, one can obtain mean Et of the corresponding distributions 
in the plots of Figs. [ri| and 12 (with the RMS values shown in the plots). 



So far we were discussing the charged decay channels of the -q meson alone. As 
for the uj — > 'k~^tt~it^ and — > vr+vr^ decays, let us note that the former is analogous 
to the 7] — > TT^^vr^vr'^ decay (even a bit easier from the viewpoint of rejection because 
the value of the charged pion deflection by angles 9 and from the vr" direction is, on 
the average, twice larger). The channel —>■ tt+tt^ does not have a neutral particle. 
In this case the probability that no signal in the HCAL will be observed after imposing 
the cut £^t|cAL-35 GeV is very small [|. 

Thus, we can conclude that to suppress contributions from the meson decays via 
charged channels one can impose the upper cuts 

• on Et deposited only in the HCAL in the isolation cone of radius -Rjsoi ~ 0.3 for 



the Barrel and for the Endcap (according to the dashed lines in Fig. [Ij). We can 
require that £'jf|"" < E'j^''. The threshold E//j;''' depends on the parent meson 
Et. We choose this cut to be equal to 0.5 — 1 GeV for Et'^ = 40 GeV, gradu- 
ally increasing it to Et^J}''=2 GeV for Et'^ = 200 GeV (see left-hand plots of 
Figs. [Tl| and 0). 

on Et deposited in the calorimeter "ECAL-i-HCAL" cells within the radius Rjgoi ~ 
0.3 beyond the 5 x 5 ECAL window. We can require that Et'^E+n'^^tE+H (with 
Et^'+H = 2 — 5 GeV/c) in accordance with the right-hand plots of Figs. 11 and 
12. These values justify our choice of the absolute isolation cut in [||] 

Therefore, we see that the two criteria introduced above may allow the charged 
hadronic decay channels of r],uj and mesons to be suppressed with a good effi- 
ciency (about 98%, depending on the exact value of Et^^^; see Figs. 11 and 12). 

We have considered the extreme case of the single meson decaying through 
charged channels. Certainly, those criteria would be much more efficient in the case of 
real pp collisions where rj, uj and mesons may have a hadronic accompaniment in 
the isolation cone around them. 

To conclude, let us add that in selection of the "7 + jet" events in real pp 
collisions these criteria can be strengthened by additionally requiring that there should 
be no of a charged particle track with Et>l GeV within the isolation cone of radius 
R]ggi = 0.7 around a photon candidate in the opposite (in 0) direction to a jet (see also 
[||], [^). Additional rejection of the charged pions that do not reveal themselves in the 
HCAL below some threshold depends mainly on the track finding efficiency for them. 

^and it is defined by the probabifity of the case when both charged pions deposit their energies in the 
ECAL alone 

^About 65 — 70% of the signal events with the direct photon satisfy this requirement (see 
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Figure 1 1 : Normalized distributions of number of the rj meson charged decay events 
over Et deposited within the radius R{'r),(/)) = 0.7, counted from the most energetic 
ECAL crystal cell, (1) in HCAL (al, bl, cl) and (2) in "ECAL+HCAL" cells beyond 
the ECAL 5x5 crystal cell window (a2, b2, c2). The 1st row (al, a2) corresponds to the 
initial Et of the 77 meson in the range 40 < -E'^nit — 60 GeV, the 2nd row to that in the 
range 60 < Et^^^-^ < 80 GeV and the 3rd to that in the range 80 < Et^^^-^ < 100 GeV. All 
distributions were obtained with the cut £^4|;^^^>35 GeV. The Barrel case (77 = 0.4). 
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Figure 12: Normalized distributions of number of the rj meson charged decay events 
over Et deposited within the radius R{r],4>) = 0.7, counted from the most energetic 
ECAL crystal cell, (1) in HCAL (al, bl, cl) and (2) in "ECAL+HCAL" cells beyond 
the ECAL 5x5 crystal cell window (a2, b2, c2). The 1st row (al, a2) corresponds to the 
initial Et of the rj meson in the range 40 < E/f^^-^ < 60 GeV, the 2nd row to that in the 
range 60 < Et^^-^ < 80 GeV and the 3rd to that in the range 80 < Et'l^^^ < 100 GeV. 
All distributions were obtained with the cut Ef^^^j^ > 35 GeV. The Endcap case 
(r? = 1.7). 
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Figure 13: Distributions of the ratios of the Et deposited in HCAL (dashed hne) and 
in "ECAL+HCAL" cells beyond the ECAL 5x5 crystal cell window (solid line) to 

the Et deposited inside this window as the function of the distance from the center of 
gravity R{r], 0), counted from the most energetic ECAL crystal cell. Left-hand and 
right-hand columns correspond to the Barrel (r} = 0.4) and Endcap (rj = 1.6) cases. 
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